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SOMMAIRE 


On  a  obtenu  dos  Equations  intdgralea  pour  les  functions  de  correlation 
dons  les  fluides  classiques  grfs.ce  A  une  methods  due  k  J.K.  Percus.  Ces  equations  sont 
appliquees  k  un  gaz  de  molecules  k  interactions  repulsives. 
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I.  rrreoTTCTTON. 

In  a  recent  paper^,  Percus  has  proposed  a  new  derivation  of  approximate 
integral  equations  for  the  pair  distribution  of  classical  fluids.  In  the  present  pap?1', 
we  shall  derive  and  discuss  some  higher  order  approximations  to  the  integral  equations 
obtained  with  Percus1  method. 


The  essential  of  this  method  is  as  follows  :  we  consider  the  distribution 
functions  in  the  presence  of  an  external  potential  iji  (r),  which  will  be  chosen  later 
on  as  the  potential  due  to  a  particle  added  to  the  system.  The  n-particle  distribution 
function  in  presence  of  the  external  potential  will  then  be  proportional  to  the  n+1 -par¬ 
ticle  distribution  function  without  the  external  potential.  On  the  other  hand,  we  choose 
some  arbitrary  functional  of  the  distribution  functions  and  of  the  external  potential 
and  we  expand  it  by  Taylor's  formula  around  ^  =  0,  the  expansion  being  made  as  an  ascen¬ 
ding  series  in  p  (r  j  (r),  where  p  (r  j  h|P  )  and  ^  (r)  respectively  are 

the  one-body  distribution  functions  in  tho  presence  and  in  the  absence  of  the  external 
potential.  With  this  method  we  can  derive  integral  equations  linking  the  various  distri¬ 
bution  functions.  The  accuracy  of  the  approximation  dopemds  on  the  number  of  terms  kept 
in  the  Taylor  expansion.  It  is  also  important  for  a  particular  problem  to  choose  in  the 
best  possible  way  the  arbitrary  functional  to  be  expanded;  that  is,  to  choose  the  func¬ 
tional  in  such  a  way  that  the  important  effects  are  taken  care  of  when  only  the  first 
terms  of  the  Taylor  expansion  are  kept;  It  remains  then  to  be  shown  that  the  contribu¬ 
tion  of  the  next  terms  of  the  expansion  is  small.  The  arbitrariness  in  the  choice  of  the 
functional  finds  its  analog  in  the  arbitrariness  in  the  choice  of  a  particular  class  of 
diagrams  which  are  summed  up  in  a  diagrammatic  approach.  It  seems  to  us  that  the  advan¬ 
tage  of  the  Percus  method  is  that,  once  the  functional  is  chosen,  the  next  approximation 
follows  automatically,  whereas  in  a  diagrammatic  approach  there  is  no  absolute  criterion 
for  selecting  which  class  of  diagrams  will  come  in  as  a  next  approximation. 


Percus  has  hown  that,  with  a  suitable  choice  of  the  functional  one  may  obtain, 

(r>\  (7,\ 

as  a  first  approximation  the  so-called  Percus-Yevi  ck '  \iT)and  Hypemetted  chain  ViT^) 


intogral  equations 
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The  first  aim  of  this  paper  is  to  derive  in  both  cases  the  intogral  equa¬ 
tions  which  rosult  when  the  next  torn  in  the  Taylor’s  expansion  is  takon  into  account 
(chapter  II).  We  call  them  for  short  FX2  and  HCN2  equations. 

Wo  then  apply  these  equations  to  particular  models  (chapter  III)  j  the 

hard  rod,  hard  square,  hard  cubo  gases,  whore  the  cluster  integrals  involving  as  much  as 

14) 

seven  points  have  been  calculated  .We  have  calculated  the  six  first  virial  coeffi¬ 
cients  for  the  PY2  equation  and  the  five  first  virial  coefficients  in  the  HCIT2  case. 

(4) 

This  may  bo  compared  with  the  exact  result  and  with  the  same  calculation  mado  with 

(5) 

the  PY  and  HNC  equation  without  correction  .  Both  the  PY2  and  the  HNC2  equations 
yield  exact  fourth  virial  coefficient.  This  is  not  so  with  the  uncorrected  equations 
except  in  the  caso  of  the  hard  rod  gas  for  which  the  PY  equation  is  exact  as  shown  in 
Appendix  I.  Whon  the  H1IC2  correction  is  taken  into  account,  the  fifth  virial  coefficient 
is  not  much  improved.  On  the  other  hand,  the  situation  is  much  more  favorable  with  the 
PY2  equation  :  the  fifth  and  sixth  virial  coefficients  remain  exact  in  tho  hard  rod  case 
and  tnko  valuos  very  near  to  tho  exact  ones  for  the  hard  square  and  hard  cube  gases. 

This  is  interesting  because  tho  hard  cubo  gas  simulates  fairly  well  the  hard  sphere  gas, 
which  itself  behaves  much  as  a  real  gas  at  high  temperatures.  As  it  is  known^that,  for 
tho  hard  sphere  gas  as  well  as  for  a  roal  gas  at  high  temperatures,  the  PY  equation  is 
a  good  first  approximation,  it  is  comforting  to  know  that  tho  next  term  in  the  Taylor 
expansion  improves  tho  PY  result  in  tho  hard  cubo  case. 


In  tho  low  tomporaturo  region,  no  conclusion  may  be  drawn  from  the  present 
study.  But,  as  it  appears  that  tho  PY  oquation  deserves  a  special  attention,  wo  huvo 
made  some  calculations  for  tho  Argon  fluid  at  low  temperatures  which  complement  tho  fow 
results  given  in  ref. (7).  It  appears  that  the  intomal  energies  and  tho  proasuros 
yielded  by  the  virial  theorem  are  practically  tho  same  as  thoso  obtained  with  the  HNC 
equation,  but  that  tho  discrepancy  betwoon  tho  pressure  calculated  from  the  virial  theo¬ 
rem  and  that  yielded  by  the  Omstein-Zomicko  relation  is  much  roducod. 

II.  DERIVATION  OF  THE  INTEGRAL  EQUATIONS. 

Wo  shall  first  recall  some  definitions  and  derive  some  properties  of  tho 
distribution  functions. 


Tho  n-porticle  distribution  function,  in  tho  presence  of  an  oxtomal  poten¬ 
tial  is  dofinod  in  the  grand  canonical  cnsombla  as  : 

. . •/ ... 
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p  (1,2  ...,n  |<J>  )  =  ^  Z  — - 

^  N  ^  n  (N-n)  I 


-/Sn 

o  dn+1  ...  dN  (i) 


whore 


is  tho  grand  partition  function  s 


^  N 

2  JL. 

N  >  0  N! 


-/®UN 


dl  ...  dN  , 


Ujj(l ,2  ...  is  the  N-particlo  interaction  energy,  which  ia  tho  sum  of  tho  interpar¬ 
ticle  potentials  V(i, j)  and  of  tho  external  potential,  ijJ  (i).  z  is  tho  fugacity. 

The  various  distribution  functions  p  (l  |<^>),  p(l,2jt|j  ),  p(l  ,2,3|i|j ) 
may  be  obtained  by  functionally  differentiating  the  grand  partition  function  with  res¬ 
pect  to  tho  external  potential  ; 


6  log 


=  e  (1  W5  ) 


/326cf'(i)6^(2) 


As^ifJk  l , 

-  J9  6  vj?  W 


:^(l,2|^)-^(lj^)p(2|ljj)+  p(l|(|>)  6(1,2) 


8  3  log  £ 

_/33^cp(i ) <5  cf)(2) 8^)  (3) 


6 2  q(i  I4>) 

/3  *  6  ip  (2)  6  qe  (3) 


=  ^  (1,2,3  jejj) 


•  •  •/ •  •  • 


I 


4. 


-  j^(  e(3  !4>)+  6(1 ,3))(p(l,2  14> )-  ^(1  (^  )  9(2|4>)) 

+  same  with  circ.  perm,  of  1,2,3^ 

-  £(l  I4>)P(2  |(j0p(3|^)+9(l|^)  6(1,2)  £>  (2,3). 


We  now  define  the  two  and  three-body  correlation  functions,  S( 1 ,2)  and  S(l,2,3),  in 


tho  absonre  of  the  external  field 

is  put  to  zero) 

:  (we  shall  in  the  f  ollowing  omit  the  1 

j)  when  it 

S(1,2)  = 

£(1,2)/(Q(1)  Q  (2)) 

(6) 

s(i ,2,3)  = 

9(1,2,3)/(9(1)9(2)(>(3)) 

(?) 

We  also  define  C ( 1 ,2)  the  direct  correlation  function  of  Omstein  and  Zcmicke  through 
the  relation  : 


-  /3  &  4*  (2) 

6  (1) 


6 (1,2) 

9(1) 


-  c(i ,2) 


(8) 


The  relation  of  a  iunctional  derivative  to  its  inverse  is  : 


/ 


d2 


&£  0) 

6 (2) 


&<±>  (2) 
6?  (3) 


&  (1,3) 


(9) 


Using  (4)  and  (8),  we  obtain  an  alternative  definition  of  C(l,2),  the  well-known 
Omstein-Zemicke  relation  : 


d3  C(1 ,3)  p(3)  (  S  (3,2)  -  1). 


S(1,2)  =  1  +  C(l ,2)  + 


J 


(10) 


5. 


Now,  if  we  take  the  external  potential  as  the  potential  due  to  a  particle 
added  to  the  system  and  put  it  at  the  point  0,  we  have  : 


poi<i>)«  iM 

P(0) 


p(l,2/>) 

P  (1,2  111))  =  -  etc... 

^  p(o) 


(11) 


Let  us  now  consider  a  quantity  f(l,2  ...  n  |i|>  )  which  is  a  function  of  the  positions 
1,2  ...  n  and  a  functional  of  if) ,  (and  thus  of  (l  |lj)))  :  as  we  have  said,  we  shall 
make  the  following  Taylor's  expansion  around  l|)  =  0  ; 


f(l ,2  ...  n  |<f>)  =  f ( 1 ,2  ...  n)  + 


( $  (n+i  |  q> )-  ?(n),;|  ^-f(1  >L.:  v.-11  1  i  l 
6<?(n+1  |  if)  ) 


lif)=C 


2! 


„  |  2 

l  dn+1  dn+2  (p(n+1  |  t())-  Q(n+1 ))( Q(r+?  (n+2))\  .£  f^’2 

J  (  A  n  fn+1/d>^(SDfn+? 


6  ^  (n+1/(p)Sp(n+2  lip 


(12) 

We  understand  immediately  how  this  procedure  is  used,  when  we  particularize  to  the  two 
functionals  considered  by  Percus  : 


a)  if  we  take 

f(l  |<J>)  -  p(l  |l|>  )  e^  , 

we  have,  using  (6),  (ll)  and  (12)  s 


(13) 


6. 


p  0  ,.Ql 


Q(0) 


.  9(1)  + 

+  jd2^(2)(S(O,2)-l)(&0.2)+^1> 


/96«p(l) 

69(2]^) 


+  1  (d2  d3O(2)^(3)(S(0,2)-0(S(0»3)-D  * 

2)  (14) 


x 


|6(i;2) 


fi&ty  (Ox 
6?(j/«|>) 


6(1.3) 


£^iiL 

69(2)4)) 


+ 


2  64>(i)__ 
6^(2 14>) 


6ct>(i)  ^ 

69 (3  | ip) 


/  >.  /3262  Ml) _ 1 

^  '^(2  14))  6p(3  I'P  )J«|>«  0 


t  two  terms  of  the  expansion,  wo  find  the  equation  : 
Tf  we  take  only  the  first  two  terms 


S(0,1)  e  "  C(0,1)  ’ 


(15) 


,  ,  4  (t0)  t0  eu^  «-  *-*  -  *->*“  C(0',)  16648  *° 

The  use  of  (15)  an  (  correlation  function  S(0,l). 

the  n  integral  oquation  fo 

«  to  obtain  the  FT2  equation. 

To  calculate  the  next  term  in  bracket  of  (U>  *  te«* 

we  must  express  the  second^derivative  to  +  .  For  this, we  use  the  relation, 

rvf  the  second  derivative  0  x 


6vp(2)64>(3) 


_  I  d4  d5  d6 


6p(il40 

6<pU) 


6p  (5 140  6p(6lij)J 
64>(2)  64>(3) 

(16) 


7. 


We  then  need,  as  may  be  seen  from  (5),  a  relation  for  the  three-body  distribution  func¬ 
tion.  This  relation  is  obtained,  as  was  done  above  for  the  two-body  distribution,  by 

e  o. *  ).*+<'>♦*<*>> 

and  keeping  the  first  two  terms  of  the  expansion  : 


£(1,2  |(|)  )  e  P  +4>(2))  =£(i,2)  + 


(17) 


r 


d3  d4  ^(3)(S(0,3)  -  1) 


6(9(1,2  |t|>)  e^(t^(l)  +t^(2)))  <£>4>(4) 


Using  (5)  and  (7),  we  obtain  the  following  linear  integral  equation  for  the  three  body 
correlation  function  which  is  expressed  in  terms  of  the  two-body  functions  S(l,2)  and 
0(1,2)  s 


S(1,2;0)  e  /3(U(1»°)  +  u(2»°)}  =  S(1,2)  +fd.3^(3)  C(0,3)(S(1,2;3)  - 


S(1,2)). 


(18) 


The  similarity  of  this  equation  with  the  PY  equation  (15),  with  (10)  taken  into  account, 
is  striking.  The  three  body  correlation  function  is  unfortunately  unsymmetrical  in  view 
of  the  special  role  played  by  the  source  particle. 


Once  we  have  the  three-body  correlation  function,  a  straightforward  manipu  - 
lation,  using  (4),  (14)  and  (16)  leads  to  the  equation  : 


S(0,1)  e^0'1)  =  S(0,1)  -  C(0,1)  +  <£  (0,1) 


• • •/ • • • 


where  <M0,l)  is  the  correction  term  calculated  as  was  said  above  s 


(19) 


8. 


$  (0,1)  =  i  jd2  d3  £(2)^(3)  C(0,2)  C(0,3) 
-  S(1,2)(1  -  e^1*2^)  -  S(l ,3)(l 


S(2,3) 


S(2,3;1) 

S(2,3) 


.  e/Svd,3)w 


d_e  /S(V(2,1)+V(3,1)) 

(20) 


Using  this  equation,  together  with  (19),  ( 1 S)  and  (10),  we  may,  implicitly,  eliminate 
(0, l),  S(2,3;l)  and  C(0,l)  and  obtain  an  equation  for  the  two  body  correlation  func¬ 
tion  which  we  call  the  PY2  equation. 


Evidently  this  approximation  process  can  still  be  carried  on  further,  but 
things  will  become  quite  complicated. 

b)  Percus  has  shown  that  if  we  take  : 


f(l  If)  =  log  (p(l|l|>  )  )  , 


(21) 


wo  obtain  as  a  first  approximation  the  HNC  equation.  Going  to  the  next  order,  we  find 
the  equation  : 


log  S(0, 1 )  c^0;l)  =  S(0,1)  -  C(0, 1 )  +  $  (°»0»  (22) 


where  the  correction  term  is  now  given  by  the  equation  : 


1  I  dif  d3£(2)  p(3)  c(0,2)  c(0,3)  s(2,3)  x 


X  .g&l?-1-)  +  1  _  log  -  S(1,2)  +  log  S(1,2)  -  3(1,3)  +  log  S(l ,3) 

S(2,3)  S(2,3) 


-  1  (  S  (0,1)  -  1  -  C(0,1))2  . 
2 


(23) 


9, 


The  equation  for  the  three-body  correlation  function  is  no  longer  linear  but  is  replaced 
by  an  equation  of  the  sane  type  as  the  HNC  equation  » 


log  o  +  V(2,0)) 

S(1,2) 


d3£>(3)  C(0,3)^Ii2i^ 


(24) 


The  combination  of  equations  (lO),  (22) ,  (23)  and  (24)  yields  the  HNC2  equation  for  the 
two  body  correlation  function. 


III.  THE  VIRIAL  EXPANSION  AND  ITS  APPLICATIONS. 

We  have  made  series  expansion  in  powers  of  the  density  for  the  FY?  and 
HNC2  equation.  The  virial  coefficients  have  been  calculated  with  the  help  of  the  compres¬ 
sibility  equation  : 

— L  s  1  -  p  (  C(r)  dr  .  (25) 

k  T  K J 


We  have  explicitly  written  down  the  relevant  Ursell  diagrams  with  their  weights  for  the 
six  first  virial  coefficients  in  the  case  of  the  PY2  equation  and  the  five  first  virial 
coefficients  in  the  case  of  the  HNC2  equation.  The  fourth  virial  coefficient  becomes 
exact  for  both  equations.  As  we  wish  to  have  some  idea  of  the  value  of  the  approximation 
for  the  next  virial  coefficients,  we  have  made  a  calculation  with  a  specific  model.  The 
hard  rod,  hard  square  and  hard  cube  gases,  for  which  the  cluster  integrals  we  need  have 
been  evaluated^  :  the  virial  coefficients  which  we  have  calculated  are  give  in  table  I 
in  the  PY2  case,  in  table  II  in  the  HNC2  case.  In  these  tables  are  also  given  the  coiv 

responding  results  of  the  exact  calculation^^  as  well  as  those  yielded  by  the  uncorrec- 
(5) 

ted  equations  , 
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Hard  rod 

Hard  squares 

Hard  cubes 

4th 

oxact 

1 

3.667 

11.333 

virial 

PY 

1 

3.778 

12.889 

coofficiont 

PY2 

1 

3.667 

11.333 

5th 

exact 

1 

3.722 

3.160 

virial 

PY 

1 

4.236 

12.434 

coefficient 

PY2 

1 

3.639 

1.701 

6th 

exact 

1 

3.025 

-18.879 

virial 

PY 

1 

4.420 

9.207 

coefficient 

PY2 

1 

3.020 

-21 .203 

Table  I.  Comparison  with  the  exact  results  of  the  fourth,  fifth  and  sixth  viiial 
coefficients  obtained  from  the  PY  and  PY2  equation  for  hard  rods,  squares  and  cubes. 

In  Table  I,  we  note  the  particular  situation  of  the  hard  rod  gas  :  in  that  case,  tho 
PY  equation  is  exact.  This  had  been  conjectured  by  Hoover  and  Poirier  and  is  demons¬ 
trated  in  Appendix  I.  The  PY2  equation  yields  also  exact  fifth  and  sixth  virial  coeffi¬ 
cients,  and  it  may  bo  supposed  that  it  is  also  exact.  The  improvement  duo  to  tho  use 
of  tho  PY2  equation  is  very  important  both  in  the  two  and  in  the  three-dimensional  cases . 
It  should  be  noted  that  in  the  calculation  of  the  sixth  virial  coefficient  we  havo  used 
in  the  correction  term  ( 20)  the  correlation  functions  calculated  with  the  lower  order 
PY  oquation,  which  is  a  consistent  procedure  in  this  approximation  schomo.  It  turns  out 
that  when  the  functions  in  front  of  tho  square  bracket  of  (20)  are  replaced  by  the  value;' 
calculated  from  the  PY2  oquation,  the  correction  is  zero  in  ono  dimension  and  amounts 
in  two  and  three  dimension  to  a  small  percentage  of  the  final  result. 
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Hard  rod 

....  _  .  . 

Hard  squares 

Hard  cubes 

4th 

exact 

- 

1 

3.667 

1U333 

virial 

HNC 

.9167 

2.944 

6.630 

coefficient 

HNC2 

1 

3.667 

11.333 

5th 

exact 

1 

3^722 

3.160 

virial 

HNC 

.8 

1 .8041 

-6.939 

coefficient 

HNC2 

1 .0835 

4.697 

—  ----- _ 

11.706 

Table  II.  Comparison  with  the  exact  results  of  the  fourth  and.  fifth  virial  coefficients 
obtained  from  the  HNC  and  HNC2  equation  for  hard  rods,  squares  and  cubes. 

The  results  for  the  HNC2  equation  are  given  in  Table  II,  Me  have  calcu¬ 
lated  only  the  fifth  virial  coefficient  and  we  may  see  that  the  improvement  brought 
about  by  the  correction  to  the  HNC  equation  is  not  very  important. 

IV.  CONCLUSIONS. 

We  have  shown  in  this  paper  that  Percus1  method  can  be  used  to  write 
down  explicitely  integral  equations  for  the  two-body  correlation  function  which  genera- 
lizeethe  PY  and  HNC  equation.  A  first  test  of  the  fitness  of  such  equations  has  been 
given  in  the  case  of  the  hard  cube  gas.  Further  information  will  be  provided  in  a 
future  paper  where  the  equations  given  in  the  present  paper  will  be  applied  to  the  case 
of  the  lattice  gas. 
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APPENDIX  I. 


The  exact  correlation  function  in  ono  dimension  is  known 


(9) 


to  be  : 


S(x)  = 


A(x  -  k)(x  -  k)k  "  1 

U-  l)k  (k-  1)! 


exp 


U  -  k) 


(A-1) 


whore 

A(x  -  k)  =  0 

if  x  4  0 

A  (x  -  k)  =  1 

if  x  s  0 

and  whore  the  length  of  the  hard  rod  has  been  taken  as  length  unit,  h  is  the  inverse 
of  the  density. 


We  define  the  following  Fourier  transform  s 


G(p) 


1 


r  +CO 

I  dx  exp  (  -ipx)  (  S(x)  -  1 ) 


J.oo 


cos  p  -  { —  1 )  p  sin  p  -  1 

(  1  -  I)2  p2/2  -  (cos  p  -  (l-  1  )  p  sin  p  -  1) 


(A-2) 


Using  (lO)  and  (A 2)  we  have  the  Fourier  transform  C(p)  of  the  direct  correlation  func¬ 
tion  : 

o(p) 

1  +  G(P) 

(A-3) 

2(cos  p  -  (  -  l)  p  sin  p  -  l) 

(4-  i)2p2 


C(p) 
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The  direct  correlation  function  is  given  by  tho  Fourier  transform  of  (A»j)  • 

1 

C(x)  =  -  - 5  (£-  x)  if  x  <  1 

U-  I)2 

(A-4) 

C(x)  =  0  if  x  >  1  . 

The  Percus -Yevick  equation  defined  by  (10)  and  ( 1 5)  is  thus  satisfied. 
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APPSNDIX  II. 


We  summarize  in  this  Appendix  the  calculations  which  have  been  made  using 
the  PY  equation  in  the  case  of  liquid  Argon  at  low  temperatures.  The  technical  compu¬ 
tation  which  has  been  used  is  the  same  as  in  ref. (7).  The  potential  is  also  the  some 

12  6 

Leonard  Jones  potential  V(r)  =  4  £  ((<s/r)  —  ( <S  /r)  ). 


The  calculation  has  been  made  for  the  following  isolated  points  : 


T  =  84°  =  0.70 

£  /k  , 

7 

d  =  1.4  gr/eni 

T  =  126°  =  1.05 

£  A  , 

d  =  1,09  gr/ cm^ 

T  =  143°  =  1.19 

&A  , 

d  =  0.87  gr/em^ 

0e84  mole/ (5  ^ • 
0.66  mole/ <5 
0.52  mole/ <3 


and  for  two  isotherms  near  the  critical  point,  at 

T  =  153°  =  1.28  e/k  and  T  =  158°  =  1 .32  g/k  . 


The  results  of  these  calculations  are  as  follows  : 

1°)  Within  the  precision  of  the  calculation,  the  pressures  calculated  from  the  virial 
theorem  and  the  internal  energies  are  the  same  as  those  reported  in  ref. (7).  This  is 
even  the  case  when,  for  the  same  density  and  temperature,  D  p/ calculated  from 
the  HNC  and  FY  equation  differ  in  sign  (this  occurs  for  a  part  of  the  158°  isotherm). 
Thus,  we  may  infer  that  the  critical  data  determined  from  the  pressures  obtained  with 
the  virial  theorem  will  be  the  same  as  in  ref. (7) (see  table  III). 

2°)  The  values  of  the  compressibility  integral  differ  from  the  samo  values  calculated 
with  the  HNC  equation.  The  discrepancy  between  the  pressures  calculated  from  the  virial 
theorem  and  from  the  compressibility  formula  is  substantially  reduced.  Wo  need  only  to 
give  one  example,  as  all  the  results  we  have  got  are  of  the  same  type  :  at  T  =  158°, 
d  =  0.60  gr/cm"1  =  365  ncle/rt  ^ ,  wo  get  for  the  compressibility  factor  P/  £  kT 
the  value  z  =  0.33  from  the  virial  theorem  (for  both  equations),  z  =  0,27  when  the 
Omstoin-Zcrnicko  formula  is  integrated  over  the  density,  in  the  FY  case  and  z  =  0,225 
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(8) 

for  the  same  in  the  HNC  case.  M.  Klein'  '  has  rightly  insisted  on  the  fact  that  there 
was  no  compelling  reason  to  prefer  the  pressures  obtained  from  the  virial  theorem  to 
those  calculated  from  the  compressibility  integral,  and  that  the  discrepancy  between 
the  two  was  somehow  a  measure  of  the  inconsistency  of  the  theory.  It  is  from  this 
point  of  view  interesting  to  compare  the  critical  data  obtained  using  the  virial 
theorem'  ; (table  III,  line  2),  and  those  obtained  from  the  compressibility  formula  : 
line  3  of  table  III  for  the  PY  case  (present  study),  and  line  4  for  the  HNC  case^. 


T 

c 

k  T 

c 

d 

c  <2 
(gr/cnr  ) 

Pc  , 
mole/o'  J 

Pc 

£ 

*ckT„ 

experiment 

151° 

1.26 

0.536 

0.324 

0.292 

virial  theorem 

150° 

1.25 

0.522 

0.315 

0.297 

f  PY 

157° 

1.31 

0.475 

0.29 

0.32 

o.z.  1 

(hnc 

167° 

1.39 

0.495 

0.30 

0.37 

Table  III.  Critical  data  : 

-  Line  1  :  experiment. 

-  Line  2  :  results  obtained  when  the  virial  theorem  is  U3od  (both  for 

the  HNC  and  PY  equations). 

-  Line  3  and  4  :  rosults  obtained  with  the  pressure  calculated  from  the 

Omstoin-Zornicke  formula  using  the  PY  and  HNC  equations. 
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